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ABSTRACT 

We test with gravitational lensing data the dark matter (DM) halos embedding the luminous baryonic com- 
ponent of galaxy clusters; our benchmark is provided by their two-stage cosmogonical development that we 
compute with its variance, and by the related 'a-profiles' we derive. The latter solve the Jeans equation for 
the self-gravitating, anisotropic DM equilibria, and yield the radial runs of the density p(r) and the velocity 
dispersion off/) in terms of the DM 'entropy' K = a 2 / ' p 2 / 3 oc r a highlighted by recent A^-body simulations; 
the former constrains the slope to the narrow range a « 1 .25 - 1 .3. These physically based a-profiles meet the 
overall requirements from gravitational lensing observations, being intrinsically flatter at the center and steeper 
in the outskirts relative to the empirical NFW formula. Specifically, we project them along the l.o.s. and 
compare with a recent extensive dataset from strong and weak lensing observations in and around the cluster 
A1689. We find an optimal fit at both small and large scales in terms of a halo constituted by an early body 
with a « 1.25 and by recent extensive outskirts, that make up an overall mass 10 I5 M Q with a concentration 
parameter c w 10 consistent with the variance we compute in the ACDM cosmogony. The resulting structure 
corresponds to a potential well shallow in the outskirts as that inferred from the X rays radiated from the hot 
electrons and baryons constituting the intracluster plasma. 

Subject headings: Dark matter — galaxies: clusters: general — galaxies: clusters: individual (A1689) — 
gravitational lensing — X-rays: galaxies: clusters 



1. INTRODUCTION 

The collisionless, cold dark matter (DM) particles that con- 
stitute the gravitationally dominant component of galaxy clus- 
ters are distributed in a 'halo' embedding the electromagneti- 
cally active baryons. The halo development under self-gravity 
from small density perturbations has been focused by several 
recent Af -body simulations, with three main outcomes. 

First, the growth is recognized (Zhao et al. 2003; Hoffman 
et al. 2007; Diemand et al. 2007) to comprise two cosmogo- 
nic stages: an early fast collapse including a few violent major 
mergers building up the halo 'body' ; a later, quasi-equilibrium 
stage where the outskirts develop from the inside-out by mi- 
nor mergers and smooth accretion (see Salvador-Sole et al. 
2007). The transition occurs at the redshift z, when a DM 
gravitational well attains its maximal depth, or the radial peak 
of the circular velocity v 2 = GM/R its maximal height (see 
Li et al. 2007). This sharp definition of halo 'formation' also 
marks the time for the early gravitational turmoil to subside. 

Second, the ensuing quasi-equilibrium structure is effec- 
tively expressed in terms of the functional K = a 2 / p 2 / 3 that 
combines the density p and the radial velocity dispersion 
a 2 , in the form of a DM 'entropy' (or rather 'adiabat'; see 
Bertschinger 1985, Taylor & Navarro 2001, Hoffman et al. 
2007, Vass et al. 2008). This mimics the behavior of a ther- 
modynamic entropy in that it increases in the halo bodies dur- 
ing the fast collapses and stays put during the subsequent quiet 
accretion. 

Third, the simple run K(r) oc r a is empirically found to 
hold in the settled halo bodies, with slopes around 1 .25 (see 
Navarro et al. 2008). This apparently universal halo feature 
allows recasting the pressure p(r)a 2 (r) oc K(r) p 5 / 3 (r) in terms 
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of the density, to balance self-gravity for the equilibrium. 

As to the latter we have used the isotropic Jeans equation to 
derive the 'a-profiles' for DM quantities like p(r) and a 2 (r), 
having ascertained that a is to lie within the narrow range 
1.25-1.3 from a state-of-the-art semianalytic study of the 
cosmogonic halo development, see Lapi & Cavaliere (2009, 
hereafter LC09). 

In this Letter we first refine these profiles to conditions 
of developing outskirts and anisotropic velocity dispersion; 
then we test them against recent data that join observations of 
strong and weak gravitational lensing (GL). 

2. HALO DEVELOPMENT AND EQUILIBRIUM 

Our framework will be provided by the standard A- 
cosmology, i.e., a flat Universe with normalized matter den- 
sity Q,m = 0.27, dark energy density fl\ = 0.73, and Hub- 
ble constant //(> = 72 km s" 1 Mpc" 1 . To bridge the matter- 
to the dark energy-dominated era across the cosmological 
crossover at z « 0.5, we shall express the redshift-time re- 
lation (1 +z) a r ? in terms of the parameter q growing from 
2/3 to 4/5. 

2.1. 77ie entropy slope from cosmogonic development 

The evolutions of the current bounding radius R, of the cir- 
cular velocity v 2 , and of the entropy K are obtained by LC09 
in terms of the halo mass M and its growth rate M from the 
simple scaling laws RocM/M 2 / 3 , v 2 ocM 2 / 3 , andK ocRM 1 / 3 . 
Whence straightforward algebra leads to express the entropy 
slope a = dlog K/dlogR as 



1 + - 



1 



(1) 



l+2e/q ' 

in terms of the inverse growth rate e = — d log (1 +z)/d logM = 
qM/Mt. 

With f«l marking the cosmogonic transition from fast to 
slow accretion as gauged on the running Hubble timescale, it 
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FIG. 1 . — Evolutionary track of a halo's mass M, circular velocity vj, entropy slope a, and concentration c for a current overall mass of 1O 15 M . In all panels: 
red lines refer to the history of an average cluster; blue ones refer to the main progenitor to illustrate the variance (see § 2.1 for details). Lines are dotted during 
the fast collapse of the halo body, and dashed during the slow accretion of the outskirts; big dots locate the cosmogonic transition, when a = 1.29 and 1.25 are 
seen to hold (bottom left panel) for the average and the main progenitor history, respectively. During the subsequent slow accretion such values are retained in 
the halo body, while they slowly decrease at the outskirts' boundary following ot(z). 



is seen from Eq. (1) that the range a « 1.25 — 1.3 will apply 
to average halos that began their slow accretion in the corre- 
sponding interval z, ~ 1.5-0.2. The range of a is narrow as 
the evolution of e(t)/q(t) is slower than for e(f) and q(t) sep- 
arately, which explains why closely similar values of a are 
found in the halo bodies from different simulations. 

Such a behavior is checked and refined in terms of the de- 
tailed evolution of e(f); its average evolutionary track is ob- 
tained from integrating forM(f) the differential equation 

[ M A 2 P 
M(M, t) = / dM' (M-M 1 ) ™2 , (2) 
Jo AM' At 

with the state-of-the-art kernel detailed in Appendix A of 
LC09. We illustrate as red lines in Fig. 1 our outcomes for 
an average cluster with current overall mass M« 10 15 M Q ; we 
plot the redshift evolution of the mass M(z), of the circular ve- 
locity v c (z), of the entropy slope a(z). Note that our approach, 
which includes ellipsoidal collapse of the body and outskirts 
growth controlled by A, renders the peaked behavior of v 2 (t) 
in remarkable agreement with the detailed simulations. For 
M ps 10 15 M Q the transition occurs at z, ~ 0.2, the entropy 
slope at Zt is around a 1.3, and the outskirts are currently 
rudimentary. 

But considerable variance arises from the stochastic nature 
of the individual growth histories. As an example of variant, 
early biased track we focus on the one associated with the 
'main progenitor' that constitutes the main branch of a merg- 
ing tree (illustrated, e.g., in Cavaliere & Menci 2007, their 
Fig. 1). Such a history obtains from the same Eq. (2), with 
the lower integration limit replaced by M/2; the results for a 
current mass M w 10 15 M Q are shown as blue lines in Fig. 1. 
Relative to the average, this history features a higher transition 
redshift z t ~ 1.5, a less massive body with M w 2 x 10 14 M Q , 
an entropy slope a w 1.25, and currently extensive outskirts. 
We find the occurrence of such biased halos relative to the av- 



erage to be 0.125 : 1 on integrating the kernel of Eq. (2) over 
the corresponding two histories. 

An imprint of the transition redshift z, is provided by the 
concentration parameter c(z), that in overall terms scales as 
[M(z) /M(z,)] 1/3 ; in fact, Zhao et al. (2003) and Wechsler et al. 
(2006) describe its increase for z < Zi with the approximation 
c(z) ~ 4(1 +Zr)/(1 + z), adopted in the bottom right panel of 
our Fig. 1. It is seen that for the average history the present 
concentration reads c ss 4, while for the main progenitor it 
takes on values c w 10. 

From overall halo development we turn now to profiles for 
the equilibrium following the transition time. 

2.2. a-profiles from Jeans equilibrium 

Physical profiles are derived from the above values of a 
inserted into the radial Jeans equation, with pressure pa 2 oc 
r a /J 5 / 3 and anisotropy described in terms of the standard Bin- 
ney (1978) parameter (3 = 1 - erg/of. In terms of the density 
slope 7 = -d log p/A log r Jeans may be recast to read 

3 v 2 6 
7= -(<*+-§) + -/?. (3) 
5 erf 5 

When supplemented with the mass definition M(< r) = 
4-n J^Ar' r' 2 p(r') entering v 2 , this constitutes an integro- 
differential equation for p(r), that by double differentiation re- 
duces to a handy 2 nd order differential equation for 7 (Austin 
et al. 2005, Dehnen & McLaughlin 2005). 

With a = const and (3 = (meaning isotropy), LC09 found 
that physical solutions, that we named 'a-profiles', exist for 
a < 35 /27 = 1 .296; the corresponding density runs steepen 
monotonically outwards and satisfy physical central and outer 
boundary conditions, respectively: a round minimum of the 
potential along with a round maximum of the pressure; a finite 
(hence definite) overall mass. In Fig. 2 we report as dashed 
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FIG. 2. — Radial runs for the a-profiles: density p, mass M, velocity dispersion o>, and circular velocity v c ; the profiles are normalized to 1 at the point 
where 7 = 70 holds (see Eq. 4), with the main body spanning the range r < 2ro. Curves are for a = 1.25 (blue) and 1.29 (red), in the isotropic (dashed) and 
anisotropic (solid) equilibria; we adopt the anisotropy given in Eq. (5) with /3(0) = -0.1 and j3' = 0.2, the latter parameter being actually irrelevant (see § 3.1). 



lines the a-profiles for various quantities: density p(r), mass 
M(< r), circular velocity v 2 c {r), and velocity dispersion aj(r). 

The behavior of the density run (top left panel of Fig. 2) is 
highlighted by the analytic expressions of the slopes 

3 3 
la = -^a, 7o = 6-3a, j b = -(l+a). (4) 

These start with the central (r — > 0) value -f a w 0.75-078, pro- 
gressively steepen to 70 w 2.25-2.1 at the point rrj that marks 
the halo main body, and steepen further into the outskirts to 
the value 3.38 -3.44. 

Monotonic behavior and physical boundary conditions for 
the a-profiles are seen (cf. LC09) to imply a maximal value 
Kcrit(a) = Vp/cr 2 ~ 2.6-2.5 in the body, at the point r p > ro 
where v 2 c (r) peaks (see Fig. 2, bottom right panel); there the 
slope reads -f p = 3(a + £ cr i t )/5 w 2.32-2.28 after Eq. (3). 

Thus the inner slopes are considerably gaffer as to yield a 
smooth central pressure, while the outer one is steeper as to 
yield a definite overall mass, compared to the empirical NFW 
formula (Navarro et al. 1997). The latter, in fact, implies 
infinite mass, and angled central pressure and potential. 

The concentration parameter for the a-profiles is given in 
detail by c = R/r-2 in terms of the radius r_2 where 7 = 2. This 
may be viewed as a measure of central condensation (small 
r-2) and/or outskirts' extension (large R). 

3 . CLUSTER PROFILES TESTED WITH GL OBSERVATIONS 

A significant observational test requires to go beyond the 
limitations to a = const and isotropy; we tackle these issues 
in turn. 

3.1. Refined a-profiles 

As to the latter, here we include on the r.h.s. of Eq. (3) the 
anisotropy term. This clearly will steepen the density run for 
positive j3 meaning radial dominance, as expected in the out- 
skirts from infalling cold matter. Tangential components must 



develop toward the center, as expected from increasing impor- 
tance of angular momentum effects (see LC09), and as sup- 
ported by numerical simulations (Austin et al. 2005; Hansen 
& Moore 2005; Dehnen & McLaughlin 2005). In detail, the 
latter suggest the effective linear approximation 

/3(r)«/3(0) + /3'[ 7 (r)- 7fl ] (5) 

with (3(0) > -0.1 and f3' « 0.2, limited to (3(r) < 0.5. We 
find (Fig. 2) the corresponding p(r) to be slightly flattened 
at the center by a weakly negative (3(0), and considerably 
steepened into the outskirts where j3(r) grows substantially 
positive. Specifically, the following simple rules apply: the 
slope (3 1 in Eq. (5) is found to drop out from the derivatives of 
the Jeans equation (Dehnen & McLaughlin 2005); the upper 
bound to a now reads a = 35/27 — 4/3(0)/27; j a is modified 
into 7 fl = 3a/5 + 6/3(0)/5 while 70 and 7^ retain their form. 

As to the former and minor issue, we include the slowly de- 
creasing run a(r) (cf. Fig. 1) enforced within outskirts devel- 
oping by slow mass accretion, as the outer scale is stretched 
out while the body stays put. Clearly, this affects little the 
inner profile of p(r) as Jeans itself (with its inner boundary 
conditions) works from the inside out, but it tilts down p(r) 
appreciably into the outskirts. 

In Fig. 2 we represent as solid lines the a-profiles refined 
as to include both the declining a(r) and the anisotropies dis- 
cussed above. 

3.2. Testing the a-profiles and their development 

Now we turn to testing these refined profiles against the 
recent, extensive GL observations of the cluster A 1689 that 
join strong and weak lensing to cover scales from 10~ 2 to 2.1 
Mpc, the latter being the virial radius R, see Broadhurst et al. 
(2005), Halkola et al. (2006), Limousin et al. (2007), Umetsu 
& Broadhurst (2008); we focus on the dataset presented by 
the latter authors and refined in Lemze et al. (2008). 

These observations concern the surface density, for which 
our benchmark is constituted by the refined a-profiles inte- 
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FIG. 3. — Surface density runs for the cluster A1689. Filled symbols represent the data by Lemze et al. (2008; see also Umetsu & Broadhurst 2008) from 
joint strong and weak GL observations. Blue lines illustrate fits to the data from our o-profiles with no prior (dot-dashed), or with priors: a > 1.25 and (3(0) = 
(dashed); a > 1.25 and (3(0) > —0.1 (solid). For comparison, the black dotted line shows the NFW fit. We report in Table 1 the values of the fitting parameters 
a, (3(0), and c and of the corresponding minimum y^/dof, 



grated over the l.o.s. at a projected distance s from the center 

rp(r) 



£0) = 2 / dr 



TABLE 1 

Results of the fits to the surface density of A 1689 



(6) 



Reference in Fig. 3 
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X 2 /dof 



In Fig. 3 we illustrate as blue lines the outcomes of fits to the 
data in terms of our refined a-profiles; for the anisotropy we 
have adopted Eq. (5), with /?(0) as the only effective parame- 
ter, see § 3.1. 

The dashed line represents the best fit for isotropic pro- 
files (J3 = 0) with the prior a > 1 .25 from the two-stage de- 
velopment (see § 2.1); the solid line refers to the best fit 
for anisotropic profiles subject to the priors a > 1.25 and 
/3(0) > -0.1 from simulations (see § 3.1); the dot-dashed line 
refers to the best fit for anisotropic profiles with no prior. For 
comparison, the black dotted line illustrates the fit with the 
NFW formula. 

In Table 1 we report the corresponding values of the fitting 
parameters a, (3(0), and c (with their 68% uncertainty), and of 
the corresponding minimum x 2 /dof. The physical a-profiles 
generally provide fits of comparable or better quality than the 
empirical NFW formula; note that an optimal fit obtains with 
no priors on a and /?, but at the cost of ignoring the informa- 
tion on the former as provided by the two-stage development, 
and on the latter as provided by the numerical simulations. 

We find that balanced fits to the surface density in A 1689 
require c w 10; lower values would cause wide overshooting 
of the outer points, while higher ones would allow approach- 
ing these (with little impact on \ 2 owing to the large uncer- 
tainties) at the cost of overshooting the precise inner points. 

A similar balancing has been found by Broadhurst et al. 
(2008) from fits with the empirical NFW formula, yield- 
ing generally high concentrations. These authors suggest 
that large values of c may be understood in terms of forma- 
tion redshifts earlier than expected from the standard ACDM 
cosmogony, as quantitatively focused by Sadeh & Rephaeli 
(2008). From the perspective of our a-profiles, concentra- 
tions c » 10 are strictly related to transition redshifts z, ~ 1 .5 
based on the state-of-the-art evolution of biased ACDM halos, 
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NOTE. — For the anisotropy we have adopted Eq. (5), where (3(0) is the 
only relevant free parameter, see § 3.1. 

see Eq. (2) and Fig. 1; on the same ground, we compute their 
occurrence to be bounded by about 13% in blind sampling. 

On the other hand, such early biased halos tend to be fa- 
vored with GL data. This is because strong GL observa- 
tions are favored by centrally flat profiles and steep outskirts 
producing conspicuously large Einstein rings (Broadhurst & 
Barkana 2008). Meanwhile, weak GL observations require 
extensive outskirts to affect numerous background galaxies. 

For A 1689 the latter data with their uncertainties are not 
yet sharply constraining the fit, but we expect convergence to- 
ward the physical profiles as the uncertainties are narrowed 
down by improved control over the redshift distribution of 
background galaxies (see Medezinski et al. 2007, Limousin 
et al. 2007, Umetsu & Broadhurst 2008). On the other hand, 
more and more clusters are being covered by GL observa- 
tions, which often find evidence of centrally flat density runs 
(Bradac et al. 2008, Sand et al. 2008, Richard et al. 2009, 
Oguri et al. 2009), consistently with our physical a-profiles. 

4. DISCUSSION AND CONCLUSIONS 

To sum up, the DM halo benchmark we test with GL ob- 
servations is comprised of a time and a space behavior, that 
we strictly link in the framework of ACDM cosmogony. As 
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to time, we find a narrow range of the DM entropy slope a as 
the outcome of a two-stage development comprising an early 
fast collapse of the body followed by a slow, inside-out growth 
of the outskirts. 

As to space, the physical a-profiles we derive feature den- 
sity runs p(r) intrinsically gaffer toward the center, and intrin- 
sically steeper toward the outskirts as to yield a definite mass, 
relative to the singular NFW rendition of early TV-body data. 
We find these runs to be stable with, or even sharpened by 
anisotropy. 

These physical a-profiles improve at both small and large 
scales the fits to the GL data, including the extensively probed 
case of A1689. Here the present analysis requires a halo 
biased toward a main progenitor lineage, with non-standard 
concentration c « 10 marking a body collapsed early at it ~ 
1.5 and late extensive outskirts; we find such halos to com- 
prise some 10% of the clusters. 

Alternative proposals to relieve the tension of the GL ob- 
servations with the profiles expected for average clusters in- 
clude density cores flattened by degenerate pressure (Naka- 
jima & Morikawa 2007), or concentrations enhanced by some 
30% owing to sharply prolate triaxialities (see Hennawi et 
al. 2007; Oguri & Blandford 2009; Corless et al. 2009) 
which on the other hand would produce steeper central slopes. 
The a-profiles and the two-stage development dispense with 
such contrasting interpretations, providing sharp and consis- 
tent shapes linked with early-biased histories. This view 
clearly invites blind sampling of more clusters in GL. 

On the other hand, as we discuss below, X-ray data will pro- 
vide an independent line of evidence concerning profiles and 
concentrations. This is based on the other major component 
of clusters, i.e., the intracluster plasma (ICP) which settles to 
its own equilibrium within the DM potential well, and emits 
strong X rays by thermal bremsstrahlung (see Sarazin 1988). 

ICP information from spectroscopy (yielding the temper- 
ature T) and X-ray brightness (yielding the squared number 
density n 2 ) are best combined in the form of the thermo- 
dynamic ICP entropy (adiabat) k(r) = ksT/n 2 ' 3 ; this modu- 
lates the ICP equilibrium within the DM potential well, and 



throughout the cluster body follows a powerlaw run k(r) oc r° 
with a < 1 . 1 . 

In fact, Lapi et al. (2005) and Cavaliere et al. (2009) com- 
pute the slope a to be expected in the cluster outskirts from ac- 
cretion of external gas shocked at about the virial R (see also 
Tozzi & Norman 2001). They find a « 2.37- 0.7 1/A0(c) in 
terms of the potential drop v^(R) A</j(c) from the turnaround 
to R; meanwhile, the ICP density follows n(r) oc r~ g with 
g sa 1.42 + 0.47/ A^(c). Values w 0.56 are seen to ap- 
ply for a-pro files with concentration c w 5, to yield a k> 1.1 
and g pa 2.2 as measured in many clusters. But in clusters with 
extensive outskirts and higher concentrations the outer poten- 
tial is shallower and Ac/)(c) smaller; when c w 10 one finds 
A0 ss 0.47, yielding an intrinsically flatter a w 0.85 (and a 
steeper g ss 2.4 in the absence of large central energy dis- 
charges). This theoretical expectation finds gratifying support 
in the flat a w 0.8 observed in A1689 by Lemze et al. (2008). 

We note that the a-profiles provide a physical benchmark 
also useful in the context of probing cold DM annihilations 
[ex p 2 (r)] or decays [oc p(r)] through diffuse 7-ray emissions 
expected from the Galaxy center (e.g., Bertone et al. 2008), 
and the positron excess recently detected by PAMELA (e.g., 
Adriani et al. 2008). We will investigate the issue elsewhere. 

Finally, since a two-stage development also applies to hot 
DM cosmogonies (though with c decreasing in cosmic time, 
see Wang & White 2008), we comment upon the case for 
cold DM on the basis of the radial run or. (r) oc p 2 ^(r)K(r) 
r a-27(r)/3 eX p eC ]; that cold DM halos will be marked by 
er 2 (r) falling down to a few hundreds km s" 1 into the outskirts, 
cf. Fig. 2 with the dynamical observations by Lemze et al. 
(2009, their Fig. 7). Such a behavior will provide evidence 
for a truly cold nature of the DM. 
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